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Little-Parks oscillation of a Mobius strip (or ring equivalently) made of a superconductor is 
studied based on Ginzburg-Landau theory. It is shown that, if the strip is wide enough, a novel 
state appears when the number of magnetic flux quanta threading the ring is close to a half odd 
integer. As a result the shape of Little-Parks oscillation of critical temperature is modified. We 
estimate the free energy of this state and give the phase diagram of the superconducting Mobius 
strip in a magnetic field. 
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Recently novel type of crystals have been grown by 
the group in Hokkaido University, 1 ' which are referred 
to as "topological matter" . These crystals are ring, disk 
or cylinder shaped with one of their crystal axis along 
the azimuthal direction. Surprisingly, even a ring of a 
quasi-one-dimensional charge-density- wave system, e.g., 
NbSe3, was created and the formation of the charge- 
density-wave was experimentally confirmed, 2 ' which in- 
dicates that these samples are good crystals and rela- 
tively free from the expected disorder originating from 
the bending of the crystal axis. 

Among these crystals we concentrate on one of the 
most subtle structure, the Mobius strip. Mobius strip 
is a representative object which has non-trivial topol- 
ogy. Therefore it is interesting to see how its topologi- 
cal nature affects the physical properties of the system. 
Experimentally, the Mobius strip of NbSe3 has been cre- 
ated. 3 ' We especially pay attention to the case where 
NbSe3 behaves as a superconductor. In this sample, the 
most conducting axis is along the azimuthal direction. 
As we see later this is actually the best configuration to 
observe non-trivial topological effects in the supercon- 
ducting state. 

The Mobius geometry was previously investigated con- 
cerning the correlated electrons in mesoscopic systems. 
Persistent current through the so-called "Mobius lad- 
der" was theoretically studied by Mila, Stafford and Cop- 
peni. 4 ' The Mobius ladder is a simplified version of 
Mobius strip which consists of only one closed elec- 
tronic channel winding twice along a circle. In case of 
a Mobius strip, however, the situation is more compli- 
cated, since there are many channels along the strip. The 
aim of this letter is to clarify how the persistent current 
(=supercurrent) through the superconducting Mobius 
strip is affected by its topological structure. 

To answer this question we study Little-Parks oscilla- 
tion, which is the oscillation of superconducting critical 
temperature as a function of the magnetic flux thread- 
ing the ring. 5,6 ' We study this phenomenon based on 
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the Ginzburg-Landau (GL) theory. Bearing the exper- 
imental situation in mind, we assume that the Mobius 
ring is formed by an array of coupled one-dimensional 
chains, although the continuum limit is taken in the end 
of the calculation. Throughout this letter, the thickness 
of the strip is disregarded and the strip is assumed to be 
two-dimensional. 
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Fig. 1. (a) Structure of the Mobius ring threaded by a magnetic 
flux <J>. (b) Chains consisting the ring. Both ends are identified. 



Our model is described by the following GL free en- 
ergy, 
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Here the ^-coordinate (0 < x < L) is taken along the 
ring. The geometry of the system is shown in fig. I . c, 
m* and — e* are the speed of light, and the mass and 
charge of a Cooper pair, respectively, v is a parame- 
ter of interchain Josephson coupling, and a and b arc 
GL parameters. The order parameter of the «-th chain 
is denoted by ipi(x). Boundary condition for the order 
parameter is given by, 



^(0) = i/>-i+i(L), ML) = ^-<+i(0). 



(2) 
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This comes from the fact that the i-th chain and the 
(— % + l)-th chain are connected, making a closed chain 
of total length 2L, which is a direct conclusion of the 
Mobius geometry (sec fig. 1). Number of chains in the 
cross section of the strip is assumed to be even (= 2K) 
for simplicity. The x-component of vector potential A x is 
set to be a constant A x = Q/L, where $ is magnetic flux 
threading the ring. In this letter we consider only the 
Aharonov-Bohm flux, neglecting the effects of magnetic 
field applied on the strip. 

We write the order parameter as ipi(x) = Xie l6i ^ x ' and 
assume the amplitude \i to be a constant. The following 
choice of 6i(x) may be natural, 



0i(x) 



irfk + 7rL,( +1 



(i<0) 



(3) 



where Zj's are real numbers such that (/j + /_j +1 )/2 is an 
integer. This makes the order parameter single- valued. 
Since the i-th chain and the (— i + l)-th chain are es- 
sentially the same chain, we assume that they have same 
amplitude and wave number: \i = X-i+i an d h = l—i+i- 
Then /j's are further limited to integers. By substituting 
these forms into Eq. (1) and performing the integration 
over < x < L, we obtain the following expression, 
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where / = <&/4>o with 0o = hc/e* . Next we minimize the 
free energy with respect to Xi's and li's. In order to min- 
imize the free energy in the bulk, namely the third term 
of the second line, we set /, — I for all i. This, however, 
gives rise to an increase of the free energy arising from 
the last term of Eq. (4), if I is an odd number. 

Here we transfer to the continuum description of the 
free energy Eq. (4) assuming that the characteristic 
length scale in the transverse direction, namely the trans- 
verse coherence length, is much larger than interchain 
spacing which we denote by a. We define the continuous 
variable by y — ai. After scaling the order parameter 
and the constants appropriately, 
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we obtain the continuum free energy as, 
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The width of the strip is given by W 

is given bya+^r (?£) (5 - /) • By setting a = 

a(0)(S- — l) and introducing zero-temperature longitudi- 



nal (|| chains) and transverse (_L chains) coherence length 
by £||(0) = h 2 /(2m*a(0)) and ^j.(0) = v/a(0), respec- 
tively, Eq. (6) can be simplified to, 
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where t is given by t + (2ir^ (0)/L) 2 (l/2 - f) 2 with t = 
T/T c and b = b' /a(0). Here T c denotes the mean field 
critical temperature without magnetic flux. In the last 
line of Eq. (6) and Eq. (7), a is kept to handle the 
diverging constant in the limit of a — ► 0. 

The solution for ^/(y) is obtained by solving the equa- 
tion, 
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Because of the assumption made just below Eq. (3), only 
the symmetric solution \I/(y) = ^(—y) is allowed. We 
should note here that the term including delta function 
is nonzero only when I is an odd number. Therefore the 
characteristics of the solutions for even-/ and odd-/ differ. 

It is easy to see that the even-/ states are ordinary 
"Little-Parks states" , where the phase of the order pa- 
rameter of each chain changes by a multiple of 27r when 
one goes around the ring once. The critical tempera- 
ture of each even-/ state is determined from t — 1 = 
and the free energy of the ordered state is obtained as 
F — — LWa(0)(t — l) 2 /(26) after a simple calculation. 

The odd-/ states are essentially different. If we go 
around the ring once the phase increases by an odd num- 
ber times 7r and, to compensate the frustration caused 
by this "fractional winding number" , a phase shift of 7r 
occurs on the line y = 0, which is actually the origin of 
the 6- function in Eq. (8). The solution for odd-/ state is 
obtained by solving the ordinary GL equation (without 
the 8- function term in Eq. (8)) in the region y > and 
y < separately, and connecting them at y = 0. The 
integration of Eq. (8) from x — — e to x — e with e being 
a positive infinitesimal number yields, 
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which gives the condition of the connection. 7 ' 

Next we estimate the critical temperature of the odd- 
/ state. Since near the critical temperature the ampli- 
tude of the order parameter is small and the coherence 
length is large, we may introduce a variational solu- 
tion ^ = i](\y\ + |) satisfying Eq. (9) with 77 being 
a variational parameter. By substituting this expres- 
sion into Eq. (7) and performing the integration over 
—W/2 < y < W/2, we obtain the variational free energy. 
The last term of Eq. (7) vanishes in the limit of a — > 0. 
The critical temperature is estimated as a temperature 
below which the solution with nonzero 77 exists. Namely, 
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1—1 + 12(£j_(0)/W) 2 = determines the critical temper- 
ature. By comparing this with the critical temperature 
of even-l state, we can determine the phase diagram near 
T c , which is schematically shown in fig. 2. As one can 
see from this figure, the odd-Z state is realized only when 
t\ > t 2 is satisfied, where t\ and t 2 are the maximum of 
the critical temperature of odd-/ state and the minimum 
of that of even-l state, respectively, as shown in fig. 2. 
Here we introduce new parameters, rj_ = £j_(0)/VF and 
r\\ = £|| (0)/L. Since t x = 1 - 12r\ and t 2 = 1 - 7r 2 rjj 
according to the present calculation, we obtain the cri- 
terion for the odd-l state as r± < irr\\ / (2v3) . 



f(°) 




Fig. 2. Critical temperature of each state labelled by I. The max- 
imum of the critical temperature of the odd-Z state is indicated 
by t\ , whereas the minimum of that of even-l state is indicated 
by t 2 . 



Next we study the low temperature region i < 1 of 
odd-l state. As a solution of Eq. (8), we adopt the fol- 
lowing formula, which satisfies \\m y _,± 00 ^(y) = const., 
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where d is determined so that Eq.(9) is satisfied. Here 
we have neglected the effect of the boundaries y = ±W/2 
assuming £_l(0) <C W. The condition for d reduces to, 
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which is fulfilled by setting d ~ a/2 because a <C £j_(0) 
is supposed. The free energy is calculated as follows. 
By substituting Eq. (10) into Eq. (7) we obtain the 
following expression, 
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Again we have taken a 
of Eq. (7) vanished. 

Next we determine the phase boundary between even- 
l and odd-l state. The free energy of even-l and odd- 
l state, normalized by LWa(0)/(2b), are given respec- 
tively by 

2 



fl 



(e) 



(1 — t) 2 + 2(1 — i)(27rri|) s 



k-f 



0{r\), 

(13) 



-(l-t) 2 + 2(l-t)(27rri|) 2 



2-' 



-^(l-t) 3 / 2 r ± +0(rf,ri 



K). 



(14) 



r(e) 



f(o) 



By equating / 2fc and f 2 k + \ {k is an integer) we obtain the 
boundary between even-l and odd-l state. For example, 
the boundary between I = and I = 1 state is given by 
/• 



V2 r± 
ott rf. 



x \f\ — t = f. From this result, we can see 

that, at T = 0, the odd-Z state can exist only when r± < 
37r 2 r 2 /(4v / 2) is satisfied, since otherwise the boundary 
between I = and I = 1 state crosses the boundary 
between I = 1 and I = 2 state at nonzero temperature, 
giving rise to direct transition between I = and I = 2 
state below the crossing temperature. 

The obtained phase diagram is given in fig. 3. Solid 
lines show the phase boundaries obtained from the 
above argument. The broken lines show guess of the 
phase boundaries in the intermediate temperature re- 
gion, where the present treatment is not applicable. As 
we see from the figure, there can be three types of be- 
haviors depending on the geometry of the system. As is 
pointed out above, the odd-l state are not at all stable 
when — ^7=ru < r i . In this case we obtain the situation 

2V3 II 2 

depicted in fig. 3 (a). If |%rjj < r±_ < ^%ry is satisfied, 
the odd-l state appears though it is not yet stable at zero 
temperature and we obtain a phase diagram like fig. 3 
(b). If r± < Tj2 r \ i s satisfied, the odd-l state survives 
until zero temperature is reached, as depicted in fig. 3 
(c). In this letter we have assumed that rn is small so that 
jj^r}, < Yjnf\\ is satisfied. In the actual Mobius crys- 
tal fabricated by Tanda and coworkers, 3 ) £y(0) is much 
larger than £j_(0) because of the sample anisotropy, and 
basically the situation (b) or (c) is rather easily realized 
even if the strip width is not so large. 

As is clarified in the above argument, the odd-Z state 
includes a 7r-phase shift line in the midst of the strip. 
This line can be viewed as a vortex line confined in the 
strip, which we call in-plane vortex line. In fig. 4 (a) 
- (d), we have schematically depicted how an in-plane 
vortex line is formed. It is interesting to see that the 
in-plane vortex line is formed from the ordinary vortex 
penetrating the strip (fig. 4 (a)), which appears in the 
intermediate state of the transition between two even- 
l states, a phenomenon usually called the "phase slip" . 
Then the state (a) must have larger free energy than the 
two even-l states before and after the phase slip. From 
this we can conclude that when the odd-l state (d) is 
more stable than even-l states, it must be more stable 
than the state (a). Although it is not clear whether or 
not the odd-l state can be, for a certain /, the most sta- 
ble configuration among all the possible configurations 
of the order parameter, we consider that the odd-/ state 
is a strong candidate because of its uniformity along the 
strip. Further investigation of this point, however, is left 
for the future studies. 

Our results are experimentally accessible using topo- 
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integer 10 ) and, therefore, it is serious for the odd-/ states. 
Because of this, the condition for the observation of odd-l 
state can be severer in actual systems. To overcome this 
point, further development of fabrication and manipula- 
tion technique of the topological matter is required. 

In summary, we have investigated the superconducting 
Mobius strip under a magnetic field. We have pointed 
out a possibility of a new state which appears when the 
number of magnetic flux quanta in the ring is close to 
a half odd integer and clarified its structure. The mag- 
netic phase diagram is obtained for systems with various 
radius and strip width. 

The authors are grateful to Prof. S. Tanda and Dr. A. 
Kanda for useful discussions. 
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Fig. 3. Phase diagram of the Mobius ring for the case of (a) 

a5s p ll < r ^ ( b ) T72 r l <r ^< 27TH and W r ^ < T72 r l 
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Odd-Z states appears in the hatched region. 
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Fig. 4. Formation process of the odd-Z state with a 7r-phase shift 
line in the midst of the Mobius strip. The strip is shown in the 
expanded form. The line AB should be identified with A'B'. The 
bold solid line indicates a vortex line. The vortex line is drawn 
also in the outside of the strip for the sake of clarity, (a) A vortex 
is penetrating the strip, (b) The vortex is stretched in the strip 
to form an in-plane vortex segment, (c) The segment is almost 
encircling the ring, (d) By annihilating the out-of-plane vortex 
segments, the in-plane segment forms a closed loop encircling 
the Mobius ring. 



logical matter. 3 ^ However, one should note that the 
Little-Parks oscillation is actually observable only for 
the samples with a diameter less than 1 /im.. 8 ^ In such 
systems, effect of thermal fluctuation is not negligible. 9 ) 
This effect is most significant when the number of the 
magnetic flux quanta in the ring is close to a half odd 



